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ABSTRACT 

We present high angular resolution PdBI images (beam of ~0.33") of the J=27— >26 line from several 
vibrational levels («7=1 and v§=\) of HC3N toward Cepheus A HW2. These images reveal the two 
main heating sources in the cluster: one centered in the disk collimating the HW2 radio jet (the 
HW2 disk), and the other associated with a hot core 0.3" northeast HW2 (the HC). This is the first 
time that vibrationally excited emission of HC3N is spatially resolved in a disk. The kinematics of 
this emission shows that the HW2 disk rotates following a Keplerian law. We derive the temperature 
profiles in the two objects from the excitation of HC3N along the HW2 disk and the HC. These profiles 
reveal that both objects are centrally heated and show temperature gradients. The inner and hotter 
regions have temperatures of 350±30K and 270±20K for the HW2 disk and the HC, respectively. In 
the cooler and outer regions, the temperature drops to 250±30K in the HW2 disk, and to 220±15K 
in the HC. The estimated luminosity of the heating source of the HW2 disk is ~2.2xlO 4 L0, and 
^3000 L for the HC. The most massive protostar in the HW2 region is the powering source of the 
HW2 radio jet. We discuss the formation of multiple systems in this cluster. The proximity of the 
HC to HW2 suggest that these sources likely form a binary system of B stars, explaining the observed 
precession of the HW2 radio jet. 

Subject headings: stars: formation — ISM: individual (Cepheus A) — ISM: molecules 



1. INTRODUCTION 



Cepheus A East, located at 700 pc (jReid et alJ l2009f> 
and w ith an IR luminosity of ~2.5xl0 L Q ( Evans et al.l 
1981). is a very active star formin g region with signposts 
of ma ssive star formation (see e.g iHughes fc Wouterlootl 
119841) . The brighte st radio continuum sou rce is the ther- 
mal radio jet HW2 (Ro driguez et al.lll994f ). which powers 
the northeast-southwest outflow seen in CO a nd HCO + 
(|Naravanan k Walker! 1 19961 : iGomez et al.lll999ft . 

Interferometric images of the molecular emission to- 
ward HW2 have revealed a very complex picture of the 
surroundings of the radio jet sugg e sting t he presence 
of a c l uster. iMartm-Pintado et all (12005ft . iPatel et all 
(12005ft . IBrogan et all (|2007ft and IComito et all (|2007l ) 
proposed that the number of sources in the HW2 system 
could be as many as five. Later, the higher- angu lar reso- 
lution observations of lJimenez-Serra et alJ (|2007ft showed 
that the molecular gas around HW2 is resolved, at least, 
into a disk centered at the radio jet (the HW2 disk), and 
an independent hot core located at ~0.4" east HW2 (the 
H C). 

IBrogan et al.l (|2007l ) reported vibrationally excited 
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emission of HC3N (hereafter, HC3N* ), previously de- 
tected bv iMartfn-Pintado et all (|2005ft . arising from an 
unresolved condensation (<1"; HW2-NE) located at the 
same position as the HC. The vibrational states of this 
molecule are excited mainly by IR radiation re-emitted 
by dust at A<50/im. The determination of the size of 
the emitting region and of the excitation temperature of 
HC3N provides a good estimate of the l u minosity of the 
heatin g object (|de Vicente et al.l 12000ft . IBrogan et al.l 
(2007) derived a temperature of 312 K for the HW2- 
NE/H C source. Assuming a siz e of ~0.6" for this ob- 
ject (|Martm-Pintado et al.|[2005t ). the expected IR lumi- 
nosity is ~2xlO 4 L0. Since this luminosity is similar to 
that measure d in H W2, and since the resolution of the 
IBrogan et al.l (|2007ft images is not high enough to dis- 
criminate between the HC and the HW2 source as the 
main heating source, the question remains whether HW2 
or the HC is the most luminous source in the HW2 clus- 
ter. 

We present high angular resolution PdBI images (beam 
of ^0.33") of the J=27— >26 rotational lines in the vj=\ 
and ^6=1 vibrational levels of HC3N, toward the Cepheus 
A HW2 region. The HW2 disk and the HC are centrally 
heated by two massive protostars. The central source of 
the HW2 disk is the most luminous object in the cluster. 

2. OBSERVATIONS AND RESULTS 
The J = 27^26 v 7 =l e and 1 f lines (£ u /fc=487K), and 
the i>6=l e and 1/ transitions (-E„/£;=883K) 7 were ob- 
served simultaneously with the PdBI in the A configura- 
tion 8 . The correlator setup provided spectral resolutions 

7 The 1)7 and v§ vibrational levels correspond to the bending 
modes of the C-C=N and C=C-C bonds. 

8 Based on observations carried out with the IRAM Plateau 
de Bure Interferometer. IRAM is supported by INSU/CNRS 
(France), MPG (Germany) and IGN (Spain) 
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Fig. 1. — Integrated intensity images of the HC3N* J=27— >26 
VT=l e (contours) and ve=lf (grey scale and thin contours) line 
emission observed toward the HC (—11 to — 9kms _1 ; upper panel) 
and the HW2 disk (from —7 to — 3 kms - 1 ; lower panel). For the 
HC, the first contour and step level are 30 (3cr) and 50mjykms — 1 
for the VT=l e map, and 10 (2c) and lOmJykms - 1 for the ve=lf 
image. For the HW2 disk, the first contour and step level are 12 
(2cr) and 12mjykms — 1 for the v?=l e emission, and 12 (2a) and 
6mjykms _1 for the vg=lf line. Negative contours correspond to 
the 3ct level. The central coordinates are [a( J2000)=22 ,l 56 m 17.98 <! 
and 5(J2000)=+ 62 o 0T49.5"l. Fille d cross shows the expected lo- 
cation of HW2 l lCuriel et al.ll2006h . The dashed lines and filled 
circles indicate the direction and positions at which the HC3N* 
excitation is calculated (Sec.|3]l. Beam size is shown at lower left 
corner. The spectra of the Vf=l e and U6=l/ lines toward the po- 
sitions of the HC and the HW2 source are reported in the inner 
panels at lower right corner. We also show the locat ion of the HC 
(filled diamond), SMA1 and SMA2 (filled squares: | Brogan et al.l 
[2007ft . HC2 and HC3 (filled tria ngles; IComito et al.feoOTI) , and R5 
(filled star; Torrefies et al. 2001) in the lower panel of this figure. 



of -80 and 160 kHz, i.e. -0.1-0.2 km s" 1 at 246 GHz. 
The synthesized beam size was 0.39" x 0.28" with a po- 
sition angle (P.A.) of 85°. We used 3C454.3 (17 Jy) 
and 3C273 (11 Jy) as bandpass calibrators; MWC349 
(1.3 Jy), as flux density calibrator; and 1928+738 (0.8 Jy) 



and 0212+735 (0.7 Jy), as phase calibrators. Calibration, 
continuum subtraction, imaging and cleaning were done 
with the GILD AS package 9 . 

The four HC3N* t>7=l and v§=\ transitions were de- 
tected toward HW2. In FigureHJ wc show the integrated 
intensity maps of the t>7=l e (contours) and ve=lf (grey 
scale and thin contours) emission from —11 to —9 kms -1 
(the HC; upper panel), and from —7 to — 3kms _1 (the 
HW2 disk; lower panel). These lines are not blended, 
and therefore were mainly used in the analysis of our 
data. Figure[T] also shows the vj=l e and vg=lf spectra 
observed toward the HC an d the HW2 source. In agree- 
ment with the SO2 images of lJimenez-Serra et al.l (|2007h . 
the HC3N* emission is resolved into two main molecular 
condensations: one centered at the HC, and the other as- 
sociated with the protostellar disk around the HW2 jet, 
the HW2 disk. For completeness, we show in FigureQ] 
(lower p anel) the location of sources SMAl, SMA2 (filled 
square s; Broaa net al.ll2007l ) . HC2, HC3 ( filled triangles; 
Comito et al.ll2007t ). and R5 (filled star; iTorrelles et ail 
2001) also reported in the region. 



2.1. The Hot Core (HC) 

The HC3N* emission from —11 to — 9kms~ 1 (upper 
panel, Figure[T]) is relatively compact and mainly arises 
from the HC. The central radial velocity of this conden- 
sation is vlsr—— lOkms -1 , and its peak emission is lo- 
cated at —0.3" nort h-east HW2. This pos i tion is similar 
to that reported bv iMartm-Pintado et al.l (|2005l ) for the 
HC, an d consistent with the location of the HW2-NE 
source (|Brogan et al.ll20"07l ). The deconvolved size of the 
v 7 =l e line emission is 0.4"x0.7" (270 AU x 480 AU), and 
for the 1)6=1/ line, 0.25"x0.4" (170 AU x 270 AU). The 
HC is likely the po wering source of the small-scale SiO 
outflow reported bv lComito et alj (|2007h toward HW2. 

2.2. The HW2 disk 

From —7 to —3 kms^ 1 , the HC3N* emission shows 
an elongated structure c entered on HW2 that rese mbles 
the SO2 disk reported bv lJimenez- Serra et all ((2007). To 
our knowledge, this is the first time that high-excitation 
HC3N* lines are spatially resolved toward a protostellar 
disk. The HW2 disk is centered at vlsr=— 5kms _1 , 
which gives a velocity difference between the HC and 
this object of — 5kms _1 . This difference has also been 
observed in the circumstellar molecu lar gas around HW2 
at larger scales (|Codella et al.ll2006l ). 

While the HC3N* v-[—l e emission arises from 
all along the disk (deconvolved size of 1.4"x0.18", 
950 AUx 120 AU), the v e = lf line is restricted to the 
inner and hotter regions closer to the HW2 source 
(size of <0.18" x 0.4", <120 x 270 AU). The orientation 
of the i>7=l e and vq=1j line emission (P.A.~ 115°) 
is roughly perpendi cular to the HW2 jet (P. A.— 46°; 
Rodrigue z et al.lll994l ). This orientation, although not 
exactly the same, is similar to that seen i n SO2, CH3CN 
and NH3 for the s ame velocity r ange (|Torrelles et al.1 
I2007D . We note that lBrogan et alj (|2007T l reported larger 
differences in the disk orientation for different molecu- 
lar tracers. However, their molecular emission images 
were obtained for a larger velocity range from ——13 to 

9 See |http://www.iram.fr/IRAMFR/ GILD AS| 
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Fig. 2. — Images of the HC3N* v?=l e emission observed toward 
HW2 for the velocity intervals from —2 to —3.7 kms -1 , —3.7 to 
-5.4 kms -1 , -5.4 to -7.1 kms" 1 , and -7.0 to -8.9 kms" 1 . The 
central radial velocities of these intervals are shown in the upper 
part of the panels. The first contour and step levels are 16 (2cr) and 
16mjykms -1 . Symbols are as in Figure[l] Beam size is shown in 
the lower left corner. 



Okms" 1 . In any case, the discrepancies in the disk orien- 
tation could be produced by either chemical or excitation 
effects. 

The integrated inten sity maps of the vr=l P lin e for 
the velocity intervals of iJimenez-Serra et al.l (|2007l ). are 
shown in Figure[2j The kinematics of this emission are 
similar to those observed from SO2 for the HW2 disk. 
The red-shifted HC3N* emission is brighter toward the 
northwest of HW2, and progressively moves to the south- 
east for blue-shifted velocities. The HC3N* emission at 
Vlsr=— 7.9 kms -1 is clearly dominated by the HC. The 
velocity difference between the north-west and the south- 
east part of the HW2 disk is —5 kms" 1 over a projected 
distance of -1000 AU. 

In Figure[3J we show the P-V diagram of the vr=l e 
emission along the HW2 disk, after smoothing the data 
to a velocity resolution of 0.76 kms -1 . We also super- 
pose the Keplerian rotation velocity curve for a star of 
18 M Q , a disk size of — 1000 AU, and an inclination an- 
gle of th e disk axis with respect to the line of sight 
of -62° (jPatel et al.l 120051 ). Although the morphology 
of the HC3N* emission is very complex, the agreement 
of the data with Keplerian rotation (thick grey line) is 
particularly good for the red-shifted part of the HW2 
disk. The blue-shifted part is co ntaminated by emis- 
sion from the HC and shocked gas (|Jimenez-Serra et al.l 
|2007|) . Therefore, the HW2 disk seems to rotate following 
a Keplerian law. The central mass of 18 M Q is consistent 
with that estimated from the source size and the HC3N 
excitation temperature for the HW2 disk ( Section [H) . 

In c ontrast with the HW2 disk scenario, iBrogan et al.l 
(2007) reported chemical differences within the molecular 
structure around HW2 at — 1"-2" scales. These authors 
proposed that this feature could be instead due to the 
superposition of two independent objects in the plane 
of the sky. Although this possibility cannot be ruled 
out, there are several observational evidences that favour 
the disk scenario: i) HC3N* (this work), SQ 2 an d NH 3 
(|Jimenez- Serra et al.1 120071 : iTorrelles et al.1 I2007T ). show 
a continuous structure whose kinematics are consistent 
with Keplerian rotation; ii) the observed HC3N* peaks 
do not seem to coincide with SMA1, HC2, HC3 or R5; for 
SMA2, this source falls —0.25" northwest the red-shifted 
HC3N* peak (we note that the uncertainty in absolute 
position is of <0.1"); hi) the dust continuum emission de- 
lineates a similar structure to that seen in HC3N*, SO2 
or NH 3 (|Torrelles et al.ll2007t ): and iv) th e spatial distri- 
bution of the H2O and CH3OH mascrs (Torrcllcs et al. 
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Fig. 3. — P-V diagram of the HC3N* 117=1,; emission along the 
HW2 disk. First contour and step level are 15 (3<r) and 10 mjy, 
respectively. Distances are c alculated with resp ect to the HW2 
source (horizontal dotted line; Curicl ct al. 2006). Vertical dotted 
lines show the central radial velocities of the HW2 disk (—5 kms" 1 ) 
and of the HC (—5 kms" 1 ). The Keplerian rotation velocity for a 
central source with mass 18 Mq, a disk size of ~1000AU, and an 
inclination angle of 62°, is also shown (thick grey line). 



119961 : ISugivama et al.ll2007t ). is consistent with the pres- 
ence of a disk with radius — 600-700 AU. Therefore, we 
propose that the observed chemical differences within 
the HW2 disk could be produced by different external 
physical conditions, or by the presence of other hot core 
sources in the vicinity of HW2. 

3. TEMPERATURE PROFILES AND HC 3 N COLUMN 
DENSITIES 

We can combine the v$=l and the vt—1 lines of HC3N 
to derive the excitation temperature of the hot gas along 
the HC and the HW2 disk (Figure[T]) by means of Boltz- 
mann diagrams. These diagrams assume LTE and op- 
tically thin emission. In case the lines were moder- 
ately optically thick, the derived excitation temperatures 
should be considered as upper limits. 

Figured] shows the radial profile of the excitation tem- 
peratures of HC3N toward the HC (filled triangles) and 
the HW2 disk (filled circles). The errors associated with 
these temperatures are <12%. For the HC, the tem- 
perature profile is centrally peaked with a maximum 
value of 270±20 K, and drops to excitation temperatures 
of <225K within the inner 0.15". This is consistent 
with the idea that the HC is internally heated by an IR 
source (see Figure 3 of Ide Vicente et all [2002). and sug- 
gests that the hot gas is highly concentrated around the 
protostar. The temperature of 27 K is l arger than that 
derived by iMartm-Pintado et all (|2005l -160 K) from 
sin gle-dish HC3N data , but lower than that calculated 
bv lBrogan et al.1 (|2007l . -312 K). In the former, dilution 
effects could account for the temperature discrepancies. 
In the latter, the excitation temperature of — 312 K is 
only an upper limit, s ince the low lying H C3N* lines are 
likely optically thick (jBrogan et al.ll2007f ). 

From the vr—l e integrated line flux (— 0.7 Jykms" 1 ), 
and assuming a temperature of —220 K, the derived 
HC 3 N column density in the HC is of -10 16 cm~ 2 . The 
partition function of the HC3N vj level has been cal- 
culated from the one in the ground vibrational state, 
but multiplied by —0.3. This corresponds to the factor 
exp[- (290X/220.FO], where 290 K is the energy of the 
fundamental rotational level in the vj state and 220 K 
is the derived excitation temperature. If we consider an 
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Fig. 4. — Radial distribution of the HC3N excitation temper- 
atures, with their errors, across the HC (filled triangles) and the 
HW2 disk (filled circles). The temperature scale for the HC is 
shown on the left, and for the HW2 disk, on the right. Arrows 
indicate the upper limits to the excitation temperatures of HC3N. 



HC3N abundance of ~ 10~ 8 (as for the Orion hot core; 
Ide Vicente et al.lfeOO^ . the estimated H 2 column density 
is ^10 24 cm -2 . Assuming a size of ~0.5", this leads to a 
circumstellar mass of ^0.6 M@ for the HC, which i s con - 
sistent with that derived bv lMartm-Pintado et al.l (|2005f ) 
from SO2. 

For the HW2 disk, the maximum excitation of HC3N 
is found at ^0.2" northwest and southeast of the ra- 
dio jet (350±40 and 365±25K, respectively), suggesting 
that the hot gas is distributed in an inner disk of ra- 
dius ~0.2" (140 AU; as estimated from FigureS]). Out- 
side this disk, the temperature falls by more than 100 K. 
From the vj—l e integrated line flux (~0.4 Jykms -1 ), 
and assuming a temperature of ^250 K, the derived col- 
u mn density of HC3N i s ~3xl0 16 cm -2 . By using Eq. 2 
of Uimenez-Serra et al.1 (|2007| ). and assuming an HC3N 
abundance 10 of ~3xl0~ 10 -10~ 9 , the estimated mass for 
the HW2 disk i s ~0.4-l M©. Th i s mass is similar to 
that o btain ed by[Pa tcl et al. (2005h. [jimenez-Serra et aTl 
(|2007t l and lTorrelles et all (j2007h . 

4. ON THE HEATING OF THE HC AND THE HW2 DISK 

In the case of radiative excitation, the HC3N gas 
temperatures are a good estimate of the dust temper- 
atures because the continuum emission of the HC and 
the HW2 disk at 20 and 45 fim is likely optically thick. 
If we assume an H2 column density of ~10 24 cm~ 2 , a 
gas-to-dust mass ratio of 100, and dust opacities of 
3360 and 1810cm 2 g~ 1 at 20 and 45 /zm for both objects 
(|Ossenkopf fc Hennind[l994T ). the derived optical depths 
are >60. This explains th e large obscu r ation seen to- 
ward HW2 at 24.5 /xm bv Ide Wit et al.l (|2009l ). Since 
gas and dust are thermally coupled, we can then es- 
timate the IR luminosity of th e central source by us - 
ing the Stefan-Boltzmann law (|de Vicente et al.l [2000) . 
From the i>7=l emission size (~0.5"), and considering 
spherical symmetry for the HC, we derive an IR lumi- 
nosity of ^3000 L© for a dust temperature of 220 K. 
Th is luminosity, which is consi stent with that obtained 
bv iMartm-Pintado et al . (2005), would corre spond to a 
ZAMS B2-type star (~10M ? : lPanagialll973h . 

If we now consider an edge-on disk geometry with ra- 
dius 475 AU and height 120 AU for the HW2 disk, we 

10 HC3N is likely being photo-dissociated in the disk. The HC3N 
abundan ce should be a factor of 10-30 l ower than that found in the 
HC (e.g. Rodriguez-Fran co et al.llT992l) . 



estimate an IR luminosity of 2.2xl0 4 L© for a dust tem- 
perature of 250 K. This is consisten t with a ZAM S B0 
star of ~18M for the HW2 source (|Panagialll973T) . 

5. DISCUSSION 

The high-angular resolution HC3N* images toward 
Cepheus A HW2 have shown that this multiple sys- 
tem contains two massive protostars. The powering 
source of the HW2 disk stands out among them with 
an IR luminosity of 2.2x10 4 Lq. The number of ioniz- 
ing photon s for a BO sta r, as inferred from the heating 
(~10 47 s~ 1 ;lPanagiall973D. is s i milar to that estimated by 
Hugh es. Cohen fe Garri ngton (1995) from the radio con- 
tinuum emission of the HW2 radio jet (^>5x 10 46 s~ x ). 
This indicates that the HW2 source is the most massive 
protostar in the region. 

The mass of the HW2 disk, as derived from HC3N, 
is relatively high (~0.4-l M©) suggesting that this ob- 
ject is at an early stage of evolution. The dy- 
namical age of the CO/H CO+ outflow (~ 5xl0 3 yr; 
iNaravanan fc Walkerl [l996) is short compared to the 
lifetime o f a photo-evaporating dis k for a ~18M Q -star 
(<10 5 vn lG""orti fc Hollenbachll2009f ). suggesting that the 
HW2 disk probably constitutes the most massive disk as - 
sociated with B stars detected so far (|Fuente et alf 2003) . 
The chemical differe nces observed within the HW2 disk 
(|Brogan et al.l l2007h . could be due to different external 
physical conditions or to other hot core sources present 
within this multiple system. 

The derived temperature profile for the HC indicates 
that this object is centrally heated by a protostar with 
an IR luminosity of ~3000L Q . External heating of the 
HC would require lumino sities of >10 5 L P) for the HW 2 
source as pointed out by IMartm-Pintado et al.l (|2005f ) . 
Shock heating could not account for the luminosity of 
the HC since the mechanic al luminosity of the outflow s 
in the region is only ~40 L© ([Narayanan fc Walkert l996) . 
Therefore, the HC hosts a massive protostar. Since this 
object has not yet ionized its surroundings, the HC is at 
an earlier stage of evolution than the HW2 source. 

Massive stars are usually found in binary systems 
(|Mason et all Il998f) . The proximity of the HC to 
HW2 (~200AU) could be interpreted as a binary sys- 
tem. Coeval formation is likely the main mecha- 
nism for the formation of l ow-mass T-Tauri binaries 
dGhez. White fc Simonl [19971 ). However, the orbits of 
these stars are preferentially ali gned with the circ umstel- 
lar disks of the primary stars (| Jensen et al.ll200~4l ). which 
contrasts with the non-coplanarity of the HC /HW2 disk 
system. Alternatively, the HC could have formed af- 
ter th e fragmentation of the HW2 disk (jKrumholz et al.l 
2009), but this would require disk-to-total-stellar mass 
ratios of ~0. 1-0.2, i.e. factors of 3 and 6 larger than that 
observed in HW2 (~0.03). 

As shown by ICunningham et aTl (|2009l ). the cap- 
ture of a massive companion (the HC) by the HW2 
disk/protostar system in an eccentric orbit could explain 
the observed precession of the HW2 jet. If we assume 
an average d orbital period of ~370 0yrs for the binary 
(Figure 7 in lCunningham et al.ll2009f ). the expected mass 
enclosed within the system would be ^25 M© for a ve- 
locity difference of 5kms _1 and an inclination angle of 
62°. This mass is similar to the sum of the masses of the 
HC and the HW2 source as derived from HC 3 N*. 
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Finally, we cannot rule out the idea that the formation 
of HW2 has triggered the formation of the HC and other 
objects in the region. Ccphcus A HW2 therefore would 
resemble the case of SgrB2, where previous episodes of 
star formation trig gered the formation of new clusters of 
massive hot cores (jde Vicente et al.ll2000[ ). 

In summary, we report the detection of HC3N* emis- 
sion toward the HC and the HW2 disk in Ccpheus A 
HW2. The HC3N* images show that these objects are 
centrally heated by massive protostars (of 18 and 10 M Q 
for the HW2 source and the HC, respectively). Since 
they appear to be very close (200 AU) , we propose that 



these objects form a binary system of massive stars. 
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